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ABSTRACT 


ADINA, a finite element program for automatic dynamic incremental nonlinear 
analysis, is used to conduct a stress analysis of the Tarawa class (LHA-1) superheater 
header. This investigation 1s sponsored bv the Naval Sea Systems Command 
(NAVSEA) as a result of cracking in the superheater header tube attachment weld. A 
2-dimensional solid axisymmetric finite element model was developed using both a 
coarse and fine mesh. Both pressure and thermaily induced stresses were studied. A 
3-dimensional solid finite element model was developed and a preliminary stress analvsis 


conducted. 
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i. INTRODUCTION 


The Tarawa (LHA-1) class amphibious assault ships have been experiencing nu- 
merous failures of their superheater header. The failures occur in the tube attachment 
weld. This is a relatively new class of ship and they are expected to be in commission 
for another twenty to thirty years. The solution of this problem is tmperative to prevent 
Beerexcessive down time experienced by these engineering plants. 

A 2 


4 


-dimensronal solid axisymmetric fintte element model was develoned using 
ADINA. A study of the effect of different thermal gradients across the thickness of the 
header was conducted. The model vrelded verv different results dependent on the 
boundary conditions imposed. A discussion ts made of the boundary conditions and the 
conclusion that a 3-dimensional model would more closely model the actual physical 


boundary conditions of this structure. 


A 3-dimensronal solid finite element model was developed. A discussion is made of 


tue boundary conditions and a study ts made of the mternal pressure loading. 


ii, DESCRIPTION OF THE PROBLEM 


NAVSEA has been conducting an extensive investigation to determine the cause for 


the weld joint cracking. The following list was provided by NAVSEA and is an analysis 


of weld joint cracking on the USS Tarawa (LHA-1): 


Linear defects were located by MT inspection in original fabrication welds and re- 
pair welds alike, including welds installed as recently as 1982. 


Predominant linear defect was the transverse Weld crack. occurring singly and in 
multipies of two or more; crack indications ranged from 1/16 inch to 3/8 inch in 
ren ciulge 

Numerous linear indications were also found across the thir L.gaments between 
welds in the A-to-D loop direction; cracks occurred singly and in multiples, and 
ranged in length from 1/16 inch to 3/4 inch. 


Predominant crack orientation was transverse to the longitudinal axis of the 
header; very few axially oriented indications were observed. 
Cracked welds in the inlet/outlet header were concentrated in the first and tome 


passes. 


No apparent concentration of defects in the welds adjacent to the division Blames 
(region of highest thermal gradient). 


The inlet/outlet header examined contained nearly 2-1/2 times as many weld joints 
with linear defects as did the intermediate header (32% of welds in no. | boner 
inlet/outlet header vs. 13% of welds in no. 2 boiler intermediate header). 


Destructive analvsis indicated the maximum crack depth to be approximately 5/32 
inch, and the crack propagation mechanism to be fatigue at elevated temperature. 


Figure 1 on page 3 and Figure 2 cn page 4 are photographs taken by NAVSEA show- 


ing magnetic particle indications of the cracking. A more detailed description of this 


problem is given in [Ref. 1: pp. 10-15]. 
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Figure 1. Magnetic particle indications of transverse weld cracks in no. 2 boiler 


intermediate header tube 





Figure 2. Magnetic particle indications of ligament cracks in no. | boiler inlet/outlet 


header tube attachment joints 


lil, DESCRIPTION OF ADINA 


ADINA is a finite element program for automatic dynamic incremental nonlinear 
analysis. The ADINA system is composed of four main programs: ADINA-IN, 
ADINA, ADINA-PLOT and ADINAT. The 1984 version was used for this thesis. The 
sequence of execution of these programs is shown in Figure 3 on page 6 which 1s taken 
mom (Ref. 2: p. 1-2]. 


A. GENERATION OF INPUT DATA 

ADINA-IN generates the input data for NDINA or ADINAT. ADINA 1s used for 
displacement and stress analysis and ADINAT performs the temperature analysis. 
ADINA-IN can be run in an interactive or batch mode. 

Coordinates can be mput in a local or global coordinate system using cartesian, cv- 
lindrical or spherical coordinates. There are several automatic mesh generation methods 
pvatiable. 

ADINA-IN has the ability to graphically check the model. All or part of the model 
can be drawn from any view desired, including element and nodal numbering. ‘Tlie 
Mists Parameter described in [Ref. 2: p. 3.25-15] gives some ability to specify which 
lines on the model will be drawn. This allows an outline of the model to be drawn and 
the ability to check for missing elements. 

The procedure for utilizing ADINA-IN to generate input data for ADINA on the 
VAX/VMS operating system is cxplained in Appendix A. The procedure for utilizing 
ADINA-IN to generate input data to solve the thermally induced stress problem ts ex- 
plained in Appendix B. The reverse Cuthill-McKee bandwidth optimization is used in 


generating the input file for ADINA. 


B. POST-PROCESSING 

ADINA-PLOT performs the post-processing of the output data from ADINA. It 
also has graphics capability, allowing a view of both the deformed and undeformed 
Siapes. the deformed scale lactorican be specificd. ADINA-PLOT can also do vector 
plotting of the principal stresses. 

Several listing options ure available. Displacements or stresses for all or part of the 
modei can be printed. Maximum values or results exceeding a specified value can also 
be listed. 
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Figure 3. Sequence of execution of ADINA prograins 


The procedure for operating ADINA-PLOT on the VAX/VMS operating system is 
explained in Appendix C. 


Sie lWOD OF SOLUTION 
ADINA uses isoparametric displacement based finite elements. The basic assump- 


Sons (rom (Ref. 3: p. 2.3-8] are: 


for the coordinates 


oe) 
ee” 


q q 
— 
sas > i . pane! | 
= Shee VS Soh z= ‘ Moy ( ‘ 
Ane — — 


for the displacements 


q q 
= a ae 7 
w= > hay — PPS w= > hy; fone! 


Peerc: 
e hf{r,s,t) = the interpolation function corresponding to node 
% 1,5.t = the isoparametric coordinates 
® q = the number of clement nodes 
¢ x,y,z, = nodal point coordinates 


® u,,v,w, = nodal point displacements 


ADINA uses the principle of virtual displacement to generate the equilibrium 


equations. The equilibrium equation as listed in [Ref. 4: 9. 124] 1s: 
AKU=R ‘Ge 


Bier: 
9 K = the element stiffness matrix 
eee — the unknown nodal point displacements 


[eee micelles vod, forces, stiriace forces, initial stresses and 


2 RK = the load vec 
ads. 


concentrated loa 


The development of the stiffness matrix for the isoparametric element is given in [Ref. 
4: pp. 203-204): 


k=| a" cBav (3.4) 


where: 


e B = the strain-displacement matrix 
ae 


ll 


the elasticity matrix 
In terms of the isoparametric coordinates: 
dV = det(/) ar ds dt (3 
ewe: 
° det(J) = the determinant of the Jacobian operator 


The stresses are calculated using the strains at the point of interest. The basic relation- 


ships used are cively im [1hGiee bp. i255 


Go 


Buy (3.6) 


r=Ce (3.7) 


@ ¢ = the element strain matrix 


© 
qq 
ll 


the element stress matrix 


The element matrices are calculated using Gauss numerical integration. A (3x3x3) tn- 


tegration order 1s utilized as specified ime eis peep 


Di AAIS NIN ETC Bie nis 


Txy = Ty, = 0 (vis the radia! direction) (3.9) 


t 


if the loading is also axisymimetric, a 2-dimensional analysis of a unit radian of the 
structure Wil provide the complete stress and strain distributions. This greatly reduces 


the computational time compared to a 3-dimensional analysis. 


E. CONTACT SURFACES 

EOIN ATaliows the specification Of contact surfaces in the development of the 
mioee!, Contact surfaces are defimed in (Ref. 3: p. 4.1-1] as surfaces that are initially in 
contact, or are anticipated to come into contact during the response solution. 

The contact surface problem is a nonlinear analysis. After cach equilibrium iter- 
ation, the segment contact force is evaluated from the nodal forces. A contactor seg- 
ment is shown in Figure 4 on page 10. The conditions of erther tension release, sliding 
@iestiching are determined for cach segment. The basic equations are given in [Ref. 3: 


p. 4.2-3]. Tension release occurs if: 
a0 (3.10) 
where: 
© = the inward normal vector to segment j as shown in Figure 4 on page 10 
@ 7) = the normal contact force for segment } 
Menmsion rclease mcans the surfaces are not in contact. 
If the surfaces are in contact the condition of sliding or sticking 1s determined by 
using the friction law. Sticking occurs if: 
eeu le ait) 
Bvere: 


— 


e Ji = the tangential contact force for segment j 


e yw = the coefficient of friction 
Sidiae occurs if: 


os 


ie. (3.12) 


The nodal forces are calculated as the consistent nodal loads which correspond to the 


tangential and normal forces of the seginent. 


@e LiGERMIAL ANALYSIS 
ADINAT uses isoparametric finite element diseretization. The basic assumptions 


aon ji.ef. 5: p. 11} are: 
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Figure 4. 
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x=) he; y=) Io 2D he eae) 
I i=l ix 
10S 
N n n 
f+ Mg=) htt Aly, t=) ity, AO= > hAO, (3.14) 
lexi i=l i= | 
where: 


® fA, = the element interpolation function 


» N = the number of nodal points 


“x, .=, — Coordinates of the nodal point 1 
e r+Ar,, = the temperature of nodal point 1 at time + Ar 
wee, «= the temperature of nodal point 1 at time t 


# AQ = the temperature increment at nodal point: 


{ 


The interpolation and substitution of these equations into the general incremental heat 


flow equilibrium cquations yield the finite clement equatrons. 


G. VERIFICATION 

Part of the installation process of a finite element code is the verification of the code. 
This is done by modeling a problem for which the results are known and comparing tlie 
Results. 

Wiiree verifications are included in this thesis. The solution of a cantilever beam 
problem is compared with classical beam theory in Appendix D. The thermal stresses 
in an annular disk are compared with classical disk theory in Appendix E. [he temper- 
ature distribution in a chimney is compared with a verified finite difference solution in 
feependix F. The last two problems were chosen because of their similarity to the 


froelen) addressed in this thesis. In each case ADINA vielded excellent results. 


I] 


IV. 2-DIMENSICNAL SOLID ANISY ViMEE TREC PiNDT Gee Eve 
MODEL 


This problem was originally modeled by Lieutenant Jon W. Kaufmann, United 


States Navy, in [Ref. 1]. The initial attempt of this thesis was to duplicate his model. 


A. INITIAL MODEL DEVELOPMENT 

NAVSEA proviced a cross section drawing of the superheater header as illustrated 
ge 13. They also provided a detail of the superheater header tube at- 
tachment weld joint in Figure 6 on page I4. The same geometry of the weld as devel- 
oped in[Ref. I: p. 1S] was utilized and is iustrategaine etre son paccal: 

A nine-node plane axisymmetric clement was used. This element is similar to the 
One Used imine 1 pens. 

A coarse and fine mesh were developed in ordemto do o study Oniconvere mec ae 
coarse mesh has 64 elements and 292 nodes. The model is shown in Figure $ on page 
16. A close-up of the weld is shown in Figure 9 on page 17, listing element numbers. 
The fine mesh has 224 clements and 1095 nodes and is illustrated tn Figure 10 on page 
18. The weld area is shown in Figure 1 on pase 19. 

The same boundary conditions stated m [Ref. |: p. 30] were apphed. The right hand 
side of the model was fixed as illustrated in Figure 12 on page 20. A frictronal surface 
Was applied between the tube and the header. The nodes on the boundary were con- 
strained together in the y direction but allowed to move with respect to each other in the 
Zz direction. This is in effect applying a contact suriace with the coefficient of friction 


equal O Zener 


B. LOAD CONDITIONS 
The same load conditions were applied as developed in [Ref. 1: pp. 38-42]. 
1. Load condition no. | - Internal pressure joad 
Load condition no. | consisted of an internal pressure of 700 psi applied at all 
nodes located on the inside of the tube and header. This corresponds to the fuil rating 


of the boiler, This loading 1s Wlasirated in ius Tr onepacer ie 
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Figure 5. Cross sectional view of the superheater header 
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Figure 6. Detail of the attachment weld joint 
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Figure 7. Geometry of the weld 
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Figure 8. 
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Figure 9. 


Close-up of the weld area, coarse mesh 
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Figure 10. Fine mesh 
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Close-up of the weld area, fine mesh 
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Fixed boundary conditions 


Figure 12. 
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Load condition no. 1 


Figure 13. 


7s). 


2. Load condition no. 2 - Longitudinal tube load 
The force per radian 1s 152 lbf. The coarse mode] uses the same consirem 
loading as [Ref. |: p. 40]. The load was applicd to the bottom element of the tubc, on 
nodes $5, 68 and 87 as indicated in the following list: 
© Fy = F’ (1/6) = 25.3 (Lb, 
¢- Ff, =F (2/5) = 101 ate 
8 F=f (1/6) = 2a) 


This is tilustrated in Figure 14 on page 23. ihe consistent loading for tle tine scsi 
applied to nodes 85, 240, 259, 238 and 87 as follows: 

@ Fs = F (1/12) = 12.66066666 (Lb, 

9 try = F (1/3) = 50.66666666 (Lb,) 

9° Pig =P GC) = 5 ocoe tee ee) 

@  Fyg = © (1/3) = 50.66666666 (Lb,) 

9 F3-= 1" (1/12) = 12.66666666 (Lb,) 


The consistent loading for the fine mesh 1s illustrated in Figure 15 on page 24. 
3. Load condition no. 3 - Thermal load 
The thermal load consists of maintaining a constant temperature of 350 degrees 
F at the interior nodes of the tube and header. A temperature gradient of [0G daame. 
F is applied to the right hand side of the model. This load condition is shown in 


Figure 16 on page: 


C. WATERTAL PROPERTIES 

The same material properties were used as listed in[Ref. I: p. 38]. Young’s Modulus 
(LE) equals 29,600,000 psi, the cocflicient of thermal expansion (a) equals .0O0Q0065 tn, 1n/ 
“F and Poisson’s rao (v) equals .3. The material of the tube and header 1s ase 
chronuum, 1% molybdenum steel (ASTM grade 22). The weld was assumied to havemiie 


same properties as the tube and header. 


De COMPARISON GPE RESULTS 

The original and deformed shapes for load conditions 1, 2 and 3 for the coarse mesh 
are shown in Figure [7 on page 26, igure {8 on page 27, and Visure 19 one 
The results for the fine: mesh are shown im Figure 20 on page 29, Fieure 21 om page 


30, and Figure 22 on page 31. The maximum stresses in the weld for the coarse and fine 
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Figure 14. Load condition no. 2, coarse mesh 
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Figure 15. Load condition no. 2, fine mesh 
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Figure 16. Load condition no. 3 
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Figure 17. Original and deformed shape of coarse mesh for load condition no. | 
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Figure 18. Original and deformed shape of coarse mesh for load condition no. 2 
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Original and deformed shape of coarse mesh for load condition no. 3 
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Figure 20. Original and deformed shape of fine mesh for load condition no. | 
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Figure 21. Original and deformed shape of fine mesh for load condition no. 2 
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Figure 22. Original and deformed shape of fine mesh for load condition no. 3 


eal 


mesh are listed in Table 1 and Table 2. For the axisymmetric problem, g,, is the hoop 


stress, o,, 1s the radial stress, and a,, is the longitudinal stress. 


Table 1. MAXIMUM STRESSES FOR THE COARSE MESH 


Stresses | Load Condition Load Condition Rada Condition 
No. | ope 


o,, (psi) -.636767E +03 S008 Or a0 - mae ae 


g,, (psi) | -.661331E+03 -.370531E+03 -.466026E +05 
-.742839E +03 -.116857E +03 | -.559437E+04 
g., (psi) | .120933E+03 | .332004E +03 .134600E +05 





Table 2. MAXIMUM STRESSES FOR THE FINE MESH 


Stresses ~ Oud Condition Load Condition | Load Condition 
No. 3 


c,,(psi) | -.649740E+03 -.534142E+03 -.136960E +06 


o,,(psi) | -.659697E+03 | .803587E+03 | -.555045E+05 







-9SOLI3E+04 
lo.psi) | .136478E+03 | .372948E ~03 | .193084E+05 


The magnitude and direction of the principal stresses for the coarse and fine mesh, 
in the area of the weld, are listed in Table 3and Table 4. A vector plot of the principal 


stresses of the coarse and fine mesh for load condition no. 3 is shown in Figure 23 on 


page 34 and Figure 24 on page 35. 


Table 3. PRINCIPAL STRESSES FOR THE COARSE MESH 





Element | PI! (psi) | P2 (psi) Actes 2) 
| -.1594E+03 -.3422E+05 | Fase 
0.1352E+04 z= -.1820E +05 | S568 
15 | 0.9153E+03 -.3779E +05 86.39 


0.3717E+04 | -.2807E +05 
0.2304E+04 | -.2084E+05 | 3618 


0.1716E+04 -.3808E+05 | $8.03 


| 0.3323E+04 -.3259E +05 


| 0.3103E+04 -.2530E +05 | $6.98 
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Table 4. PRINCIPAL STRESSES FOR THE FINE MESH. 
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The maximum hoop stress listed in [Ref. 1: pp. 39-40] for load condition no. 1 is .69 
(ksi) and the maximum hoop stress for load condition no. 2 is .96 (ksi). Those results 
compare favorably with the results shown here from ADINA. [Ref. 1: p. 42] indicates 
the maximum hoop stress for load condition no. 3 is 19.0 (ksi). This value is much 
smaller than the result of 137 (ksi) obtained with ADINA. In considering the verifica- 
tion of ADINA as explained in chapter three, it is believed that the results obtained here 


are correct. The results confirm that the thermal stresses are the predominant stresses. 
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Picurer23: 


Vector plot of the principal stresses of the coarse mesh for load condition 
now 
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Figure 24. Vector plot of the principal stresses of the fine mesh for load condition 
no. 3 
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The accuracy of the results of the problem that has been modeled are very good. 
There is very little difference between the results of the coarse and fine mesh which 
means convergence has been achieved. There is no significant difference in the calcu- 
lated stresses at the boundaries of adjacent elements. 

The next step was to use the contact surface option of ADINA for the boundary 
between the tube and the header. This allows an indication of the coefficient of friction 
and a determination of its effect on the solution. The maximum stresses are shown in 
Table 5 and Table 6 for the coarse and fine meshifera cocificient of Mictionmeieae: 


These results compare vere well wilt the previous pesiiis: 


Toavies Ss. MAXIMUM SURESSES FOR. THE VGGARSE MESH. COR 
TACT SURFACEIOV TION 






Stresses Load Condition Eoad Condition Load Condition 
Nigel Non Sore) 


lac. (psi) | -.646923E +03 -.S]4228E +03 -.] 38265 +06 
| a, (psi) -.GSTOI4IE +03 362559E +05 -465231E +05 
a.(psi) | ~747963E+03 -112500E +03 | -.S58613E+04 
| a.. (psi) boo ee aan REN 5b ie FE Ag io Sus 


Table 6. MAXIMUM STRESSES FOR THE FINE MESH, CONTACT 
SURFACE OPTION 















Stresses 






Eoad € onetien 
KO > 


si) -.646745E +03 537235E+03 7053E +06 
on ae -.657250E +03 495395E +03 
y_Pnrstoare 03 [1399096 03 


GA O51) sos 0 Us S00 704r 40s OS URE cette 


Load Condition 
INO 


Load Condition 
Non! 








E.. FORTHER MODEL DEVELOP ie 

A detailed study of the problem and the results obtained in the previous section tn- 
dicated some changes were needed tn the model. 

It was discovered that the model as developed was not a repeating section wee 
distance between the center of the tubes ts two inches. The right hand side of tlre mame 
is two inches from the center of the tube. A new model was developed with the right 


hand side berng one inch from the center of the tube. The coarse mesh consists oie 
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elements and 211 nodes. The fine mesh has 172 elements and 763 nodes. The models 
are shown in Figure 25 on page 38 and Figure 26 on page 39. 

Figure 19 on page 28 and Figure 22 on page 31 prompted a change in boundary 
conditions. The deformed shapes showed that the header was not allowed to expand 
uniformly in the axial direction. The bottom node of the right hand side of the header 
was fixed and the rest of the nodes on the right hand side were allowed to roll in the z 
direction. These boundary conditions are shown in Figure 27 on page 40. 

The original and deformed shapes for load conditions |, 2 and 3 for the coarse mesh 
ape shown in Figure 28 on page 41, Figure 29 on page 42 and Figure 30 on page 43. 
Figure 30 on page 43 shows that the header is now expanding uniformly in the axial 
direction. The maximum stresses are indicated in Table 7 and Table 8 for the coarse 


and fine mesh. The ADINA-IN input file for these models 1s listed in Appendix G. 


Mee 7. NiANTMUN STRESSES FOR THE COARSE NIESHE, RE- 
PEA TING SECTION 
Stresses Bead @endition | Load Conditions Load Condition 
iow! land 2 Nor} 


me (151) - 420381E+03 - J|60508E +03 = eg te OD 


-.597283E +03 TINIE +03 -.373304E +05 
-.§25134E +03 -.112852E +04 -454326F +04 


eel eet) AOR GO 2 05 





Table 8. MAXIMUM STRESSES FOR THE FINE MESH, REPEATING 












SECTION 
Load Conditions oed Gonaiion 
Ron! lana INOwS 
G,, (psi) -A7ZTSIE+03 - SIDSITSE + 03 - LOIJO3ZE +06 
- S40SS2E+03 ee Solely 
1 G.. (Psi) -.922094E +03 -.126150E +04 


o,, (DS1) Ploy 2 lea) O02645E +04 AS1684E + 05 


Pe YHERMAL GRADIENTS. 
the previous change in the model made the thermal gradient in the area of the weld 


very flat. It was desired to determine what effect different thermal gradients on the right 
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Figure 25. Coarse mesh, repeating section 
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Fine mesh, repeating section 


Figure 26. 
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Figure 27. 
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Figure 28. Original and deformed shape of coarse mesh for load condition no.1, re- 


peating section 
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Figure 29. Original and deformed shape of coarse mesh for load condition no. 2 , 


repeating section 
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Figure 30. Original and deformed shape of coarse mesh for load condition no. 3 , 


repeating section 
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hand side of the model would have on the results. The maximum stresses for a thermal 


sradient of 100, 10, 2 and O'decnees F are Shown mr Pable 7) 


Table 9. MANIMUM STRESSES FOR DIFFERENT THERMAL GRADI- 


ENDS 
Thermal 
Gradient o,, (psi) c,, (psi) a., (psi) By (eR, 


100°F Se] elle) Se ih : -454326E+04 | .104692E +05 


Ole 99S) Oils a ne -4A54017E +04 | .L0d558E+05 


| 2° -.998993E +05 | -.373106E +05 | -.453990E + O41; 2910 ome 
lor -.999D1ISE+05 | -.373102E+05 | -.453983E+04 | .104321E +05 


The results were surprising. The size of the gradient had ven hittle cilecuia amare 








magnitude of the stresses. The stresses remained very high even for a zero thermal gra- 
dient. A study of these results indicated the need for different boundary conditions. 
Restricting the nght hand side of the model to move in the z direction only was not al- 
lowing the header to expand outwardly in the radial direction. The header was being 
forced to expand into the center of the tube which was creatiig the high Stressesumeeme 
area of the weld. 

The boundary condition chosen 1s shown in Figure 31 on page 45. The bottom left 
hand node of the header at the boundary of the tube and the header was fixed” Titeiess 
of the bottom nodes of the header were allowed to roll in the yv direction. [he comma 
surface option of ADINA could not be used for this solution. The nodes on the 
boundary Were constrained together in the v direction and allowed to move with respect 
to one another in the z direction, modeling a contact surface with zero coefficient of 
fricuon. This was with the exception of the bottom Ic{t hand node of the header. 

An iterative procedure was then followed to determine the solution of the proSigim 
After each run of the problem, the bottom left hand node of the header was sivenme 
permanent displacement to the right equal to the displacement of the corresponding 
node of the tube at the boundary. This was continued until no overiap occurred.” ie 
Original and deformed shapes of the final iteration for load condition no. 3 are shown in 


Figure 32 on page 46. As can be seen, the model is now expanding uniformly in the 


o4 


Figure 31. 
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Figure 32. Original and deformed shape of the fine mesh for load condition no. 3, 


expanding radially and axially 
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axial direction and outwardly in the radial direction. The maximum stresses in the weld 


area of the coarse mesh for load condition no. 3 are shown in Table 10. 


Table 10. MAXIMUM STRESSES FOR 
THE COARSE MESH, EXPAND- 
ING RADIALLY 


Load Condition No. 3 
o., (psi) 460123E+03 
-.932915E +02 


se E03 
moose 1Oe +02 













G,, (psi) 


These results are much lower than those obtained for the previous model and it is 
ipeeved that these results are closer to the actual stresses. However, it is also believed 
that this model is not a true repeating section in that the expansion of the neighboring 
tubes and the rest of the header will have an effect on this section, probably raising the 
stresses in the weld area. The need for a 3-dimensional model was indicated. 

A 3-dimensional model was also desired in order to apply non-symmetric loading to 
the model. The major difference between this header and other headers with a similar 
cross section is that this header is not sectionalized. It has one inlet/outlet header and 
One intermediate header. This results in a larger temperature gradient from the inlet to 
the outlet side. A 3-dimensional model would allow this type of non-symmetric loading 
to be studied. 
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V. 3-DIMENSIONAL SOLID FINITE ELEMENT MODEL 


A. MODEL DEVELOPMENT 

The repeating section chosen to be modeled was one tube row. Only half of the row 
was modeled due to symmetry. The welds were assumed to have the same material 
properties as the tube and header and no attempt was made to substructure the weld 
area. 

The element chosen to mode! the structure was the 20-node tsoparametric sohid ele- 
ment. This element is shown in Figure 33 on page 49. 

Figure 34 on page 50 shows the modeling of one tube. Each tube consists of sixty 
elements. Figure 35 on page 5] shows now the header was connected to tiesto ema. 
contact surface was placed between the tube and the header. The entire model is shown 
in Figure 36 on nage 52. Figure 37 on page 33 is another view of the model Simgaaine 
the location of the origin of the axis. The mode! consists of 288 elementsiacuee es 
nodes. Figure 38 on page 54, Figure 39 on page 55 and Figure 40 on page 56 show a 
view of the model looking from the x, y and z directions. 

The hidden line computer code developed by David R. Hedglev, Jr. of the Ames 
Research Center was utilized to check the model. Diflerent views are sh@mgaman 
Figure 41 on page 57 and Figure 42 on page S58. 

The following boundary conditions were imposed. Every node in the y=0 plane 
was constrained to roll in that plane. In addition the node at the origin was fixed and 
the nodes above the origin along the positive X axis Were constrained to roll only in the 


xX direction. 


B. INTERNAL PRESSURE LOAD 

The internal pressure load of 700 psi was applied to the inside nodes of the tube and 
header. [t was then discovered that ADINA does not optinuze the band-width of con- 
tact surfaces. The bandwidth of this model is approximately 2700 and the problem was 
too large to be solved on the VAN; VMS operating svstem. The ADINA-IN input file 
for the contact surfaces is listed in Appendix H for reference. 

A contact surface was then modeled by constraining the corresponding nodes of the 
tube and header to move together in the y and z direction, but ailowing them to move 
relative to one another in the x direction. This assumes a coefficient of friction of zero 


and reduced the bandwidth to approximately 1000. 
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Figure 33. 20-node isoparametric solid element 
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Figure 35. Attachment of header to a tube 
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Figure 36. 3-dimensional mode! 
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3-dimensional model showing the location of the origin 
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View of 3-dimensional model from the x direction 


Figure 38. 
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Figure 39. View of 3-dimensional model from the y direction 
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View of 3-dimensional model from the z direction 


Figure 40. 


56 





Figure 41. View of model using hidden line removal 
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Figure 42. Another view of the model using hidden line removal 
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It then became necessary to enforce symmetry on either side of the z=0 plane. 
Every node in this plane was constrained to roll in the plane. It is understood that this 
enforced symmetry is more costly than modeling only half of the structure. The entire 
section was modeled to allow a study of unsymmetric temperature loading. The 
ADINA-IN input file is listed in Appendix I. 

Figure 43 on page 60 shows the original and deformed shapes for this loading con- 


dition. Table 11 lists the maximum stresses in the area of the weld. 


Table 11. NIAATNIUM STRESSES 
3-DIMENSIONAL MODEL 


| .184639E +04 
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Figure 43. Original and deformed shape of the 3-dimensional model, internal pres- 


sure load 
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Vi. CONCLUSIONS 


A. DISCUSSION OF RESULTS 

The maximum stresses obtained for the axisymmetric model were between 137 (ks1) 
and .46 (ksi), dependent upon the boundary conditions imposed. The axisymmetric 
model does not adequately represent the physical boundary conditions and the actual 
Sitess tS Somewhere between these two figures. 

A 3-dimensional model has been developed which allows an accurate representation 
of the boundary conditions. The maximum stress for the internal pressure load in the 
area of the weld was 1.85 (ksi). This model also allows non-symmetric loading to be 


applied. 


B. OPPORTUNITIES FOR FURTHER RESEARCH 

The next step in solving this problem is to obtain results for the thermal loading. 
The initial loading should be static but the transient loading should also be studied. In 
Sewas iis the transient problem, a consideration sliould be made of the temperature de- 
pendence of the material properties of the structure. Another area for research would 
be the consideration of creep. 

The major attempt of this report was to develop a 3-dimensional model which could 
be used to investigate the farlures that have occurred tn the superheater header attach- 
ment welds. This model can also be used to investigate any needed design changes. It 


is hoped that this mode! will be used for further investigation of this problem. 
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APPENDIX A. ADINA-IN PROCEDURES 


These procedures are for the VAX/VMS operating system utilizing the plot-10 
graphics library. 

To operate ADINA-IN in the batch mode, a data file which defines the problem 
must be created. Assign this data file to a logical unit number and run ADINA-IN. 

Now create the database. All of the input to ADINA-IN its stored on this danibase: 
The Command is DATAB CREATE. Next read the data file by using the command: 
READ (logical unit number). 

Tiie model can now be plotted using the frame, mesh and plot comimands as ex- 
plained in [Ref. 2: pp. 3.25-1 to 5.25-20]. The CONTROL command must be Semiemmae 
correct parameters to use graphics. For the VAX, VMS operating sysicm utilizing the 
plot-10 graphics library, set baudrate=4$00, svstem=2, device=2 and option= 1024. 
Changes can be made to the model by exiting and changing the Jata file er by typing in 
changes interactively. All chamees will be si casamine diataaase. 

When the problem 1s properly defined, create the ADINA input file using the com- 
mand ADINA. Exit ADINA-IN by tvping in END. 

ADINA-IN will have created two files. FOROOQI.DAT is the database, [litera 
base can best be saved by renaming it to a umique file name. To make changes assign 
the database to logical unit FOROOi and run ADINA-IN. Use the command DATA- 
BASE OPEN and then make any necessary ciances: 

The second file created bv ADINA-IN is FORO02.DAT and this is the input file for 
ADINA. This file should be renamed to a unrque filename and assigned to logical unit 
number FOROOS. Assign an output file to logical unit number FOROO6. Ensure there 
are no FOR*.* files in the directorv as this can interfere with the operation of (WD ia 


Now run ADINA. The requested results will be printed mn the assigned output file. 


APPENDIX B. ADINA-IN PROCEDURES FOR THE THERMAL STRESS 
PROSE 


ADINA-IN is also used for creating the input file for ADINAT. This particular 
version of ADINA-IN will not create the entire file which necessitates some minor ed- 
iting of the ADINAT input file before ADINAT can be run. 

The temperature loading in this thesis was done by using the LOADS TEMPER- 
eemeie command in the ADINA-IN input file. After the problem is defined in 
ADINA-IN, use the ADINAT command to generate the ADINAT input file. 

ADINA-IN again creates two files. FOROOL.DAT is the database and 
FOROO2.DAT is the ADINAT input file. Rename these files. 

The ADINAT input file must now be edited. Search for the line which has the let- 
fess SKARKKWK. Thus line must be deleted and the specific heat entered. The line just 
@eove this line must also be deleted and the thermal conductivity cittered. Exit tlic 
EOL NAT input file. 

Assign the ADINAT input file to logical unit number FOROUS and assign an output 
fe to logical unit number FOROO6. Run ADINAT. ADINAT creates several files, 
most of which are used for the operation of the code. The final nodal point temperatures 
and heat fluxes are printed in the output file. The nodal temperatures are also stored in 
FORO56.DAT. This file must be renamed to a unique filename. 

ADINA-IN must now be run again to create the input file for ADINA to solve for 
the thermal stresses. Assign the nodal temperature file to logical unit number FOROSO. 

ihe ADINA-IN mput file must be changed. Delete the LOADS TEMPERATURE 
command and indicate on the MASTER command line that ITP56=1. Run 
aati and read thie changed input file. Ensure all FOR*.* files are deleted before 
doing this. Then use the ADINA command to create the ADINA input file. ADINA 


can now be run as explained in the last paragraph of Appendix A. 
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APPENDIX C. ADINA-PLOT PROCEDURES 


ADINA created several files besides the assigned output file. Most of these files are 
used for the operation of the code. File FOROGODAT is the porthole tiem 
ADINA-PLOT, 

Before running ADINA-PLOT, ensure all FOROOL.DAT files are deleted and that 
logical unit number FOROOI! 1s deassigned. Run ADINA-PLOT, First create the da- 
tabase using the DATABASE CREATE command. All of the information in the 
porthole file 1s now loaded. 

Ensure the control information ts correct. Tlie model can now be plotted and the 
results listed as explained in [Ref. 6]. 

Exit ADINA-PLOT with the command END. ADINA-PLOT will have created a 
database with hlename FOROOI.DAT. Rename this file. To look at the resulisia aa 
assign the database to logical unit number FOROO! and use the DA TABNSE Reis 


conunand when running ADINA-PLOT, 
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APPENDIX D. VERIFICATION OF ADINA - CANTILEVER BEAMI 
PROBLEMNI 


The cantilever beam problem modeled is shown in Figure 44 on page 66. A moment 
of 360 k-in is applied to the end of a cantilever beam which 1s 48 inches long and 12 
inches wide. 


Wine classical Equations as listed in (Ref. 7; pp. 215 and 7/37] are: 
a= \y/I (4.1) 
5, = M,L*/2EI AO) 


The finite element model showing both the original and deformed shapes is shown in 
Figure 45 on page 67. A nine node plane stress element was used. Table 12 shows a 
comparison of the results for several nodes. As can be seen. ADINA yields exact results 
for this problem. The following is the ADINA-IN input file: 


Table 12. COMPARISON OF RESULTS 


a ) 
stress atnode Pips «dS —SSCSC~S IS 










DATAB CREATE 

Pea CANTILEVER BEAM 

ras 

MASTER REACTIONS=YES IDOF=100111 
PRINTOUT MAX 

PVaAn SES TYPE=STATIC 


* 

COORDINATES / ENTRIES NODE KX Y 4 
A C0 6 
5 0 48 0 
AUS: QO 48 12 
Jak Cage le 

a 
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Figure 44. Cantilever beam with concentrated moment 
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Figure 45. Original and deformed shapes for the cantilever beam 
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MATERIAL 1 ELASTIC E=3E4 NU=0.0 D=l 
* 


EGROUP 1 PLANE SUBTYPE=STRESSZ RESULIS=STRESSES sali —c 
STRESSTABLE N=1 Pl=1 P2=2 P3=3 P4=4 P5=5 Po=6 P7=7 P8=8 PI=9 
EDATA / ENTRIES EL THICK 

De get 

21 
GSURFACE 1 5 15 22) ERls2 EE2-2) Nopes=2 
as 


BOUNDARIES 111111 TYPE=NODES / 1 6 11 
* 

LOADS CONCENTRATED 

552650 

BO: 27S 50 

* 


CONTROL B=4800 SYSTEM=2 DEVICE=2 OPTION=1024 
CONTROL PLOTE=W PLOTS=N 
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APPENDIX E. VERIFICATION OF ADINA - ANNULAR DISK 


PROBLEM 


The classical equations for the thermal stresses are listed in [Ref. 8: p. 259]: 


r b 
Gece -uP| Trdr + (r? — a’\r7(b? - °)| Trdr] (1) 
a a 
ee 2 Ce Sa el 
og=al[—T+ ifr) Trdrt+(r +a )[r' (bo —a’)| Trdr] (E.2) 
a a 
where: 
®ha = inner radius, .63 inch 
© b = outer radius, 3.0 inch 
Pe eniperature disttioution, 399.569212 - $6.5436r + 12.160239867r° 
2 E = 29600000 psi 
¢ o¢ = 0000065 in/inj°F 


The finite element model showing the original and deformed shapes is shown in 


Figure 46 on page 70. A nine node plane axisymmetric element was used. Table 13 


shows a comparison of the results. The comparison is very good. Most of the difference 


can be attributed to the fact that the results for the disk theory were computed at the 


nodes where as the results for ADINA are calculated at the Gaussian integration points. 
The following is the ADINA-IN input file: 


Table 13. 


r (inch) 


08) 
mSl5 
feleol} 
j 3.0 


COMPAKISON OF NES UEIS 


Stresses from Disk Theory 
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Stresses from ADINA 
(pst) 


= - 120961E+05 
eee Ole Ue 
1L69057E + 04 

a, = .1LI644E+04 





DEFORHED ct —4 ).00505 


ORIGINAL wv ou D.£32 
TIME 1.000 


DISC 





Figure 46. Original and deformed shapes for the disk 
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PATAB CREATE 

HEAD °“DISC’ 

CONTROL B=4800 SYSTEM=2 DEVICE=2 OPTION=1024 
CONTROL PLOTE=WAIT PLOTS=NO 

ral 


COORDINATES / ENTRIES NODE X Y Z 


1 0 63 0 
Zz 0 63 S 
3 0 BO -<0 
x 0 S40 3 


* 

MATE REAL 1 THERMO-ELASTIC 

*™2U00 29600000 .3 .0000065 
ZV00 29600000 .3 .0000065 

* 


EGROUP 1 PLANE SU=AXISYMM 
GSURFACE 4 2 1 3 EL1=4 EL2=1 NO=9 
* 


MASTER REACTIONS=YES IDOF=100111 
PRINTOUT MAX 

ANALYSIS TYPE=STATIC 

x 


BOUNDARIES 1OLE11 TYPE =NODES/12 STEP 1 TO 20 
* 


LOADS TEMP TR=0 


oe 520 
Pa 350 
Zo 350 
Bee 250 
4 250 
2 250 


* 


cal 


APPENDIX F. VERIFICATION OF ADINA - CHIMNEY PROBLEM 


The chimney problem modeled is shown in Figure 47 on page 73. The inside surface 
is at a uniform temperature of 650 degrees C and the outside surface is at a temperature 
of 150 degrees C. 

The finite element model is shown in Figure 48 on page 74. A nine node plane 
Stress element was used. 

A finite difference temperature distribution of this problem was obtained from Pro- 
fessor G. Cantin of the Naval Postgraduate School. The results had convergequioumn. 
exact solution. Table 14 shows a comparison of the results for several nodes. The 


comparison is very good. The following is the ADINA-IN input “le: 


Table 14. RESULTS FOR CHIMNEY PROGiEEM 


Temperatures from Finite {| Temperatures form 
Ditierence ©) ADINA (C) 


| 3.6933026E +02 3.68497E +02 


3.4661577E +62 3.39891E +02 


eae > 5149076E = 02 7 52863E = 02 
BS — 3.480426 +02 


| 4.01579E +02 





DATAB CREATE 

HEAD “CHIMNEY’ 

CONTROL B=4800 SYSTEM=2 DEVICE=2 OPTION=1024 
CONTROL PLOTE=WAIT PLOTS=NO 

a 


COORDINATES/ENTRIES NODE X Y Z 
aE Os 20> el 
Z Ola OR 
S Oy 1 30) Miao 
4 QO 2.0 fio 
5 0 0 ea) 
6 0 O 0 
7 Oe 220 no 


* 

MATERIAL 1 THERMO-ELAS TIC 
-2000 29600 .3 .0000065 
2000 29600 .3 .0000065 


(2 





Figure 47. Chimney problem 
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Figure 48. 





Finite element model of the chimney 
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EGROUP 1 PLANE SU=STRESS2 
EDATA 

14 

ral 

3} 1 

GSURFACE 1 2 3 4 ELl=1 EL2=1 NO=9 
GSURFACE 4 5 6 7 ELl=2 EL2=1 NO=9 
* 

MASTER REACTIONS=YES IDOF=100111 
PRINTOUT MAX 

ANALYSIS TYPE=STATIC 

* 


LOADS TEMPERATURE 


Ll Iksyo) 
2 150 
4 S50 
14 150 
7 150 
19 aye, 
20 ILS 
21 50 
2 650 
ot 650 
3 650 
is 650 
5 650 


He 


PROBLEM 


DATAB CREATE 


*xINPUT FOR THE COARSE MESH, 2-DIMENSIONAL AXISYMMETRIC MODEL 


HEAD ‘REPEAT HEADER’ 


CONTROL B=4800 SYSTEM=2 DEVICE=2 OPTION=1024 


CONTROL PLOTE=WAIT PLOTS=NO 
*COORDINATES ARE DEFINED 


COORDINATES / ENTRIES NODE X a6 
at 0 53 
3 0 093 
oy) 0 63 
SI, 0 215 
60 0 oS 
61 0 -7D 
70 0 762 
Tek 0 Ee: 
72 0 - 75001 
85 0 263 
a7 0 iss) 
100 0 2255 
200 0 saz 
152 0 -3358 
204 0 ~9a2 2366 
PAaUpl 0 201 
210 0 - 84513864 
207 0 son 
209 0 -813 
126 0 -9o0 
178 0 3.938 
208 0 2610 
206 0 70 
50 0 1.0 
6 0 1.0 
7D 0 0 


AMATERIAL PROPERTIES ARE DEFINED 
MATERIAL 1 THERMO-ELASTIC 

-2000 29600000 .3 .0000065 

2000 29600000 .3 .0000065 
ATHE ELEMENT IS DEFINED 
EGROUP 1 PLANE SU=AXISYMM 
*THE MESH FOR THE TUBE IS GENERATED 
GSURFACE 3 1 37 39 EL1=1 EL2=3 NO=9 NC=B 
GSURFACE 39 37 60 61 EL1=1 EL2=2 NO=9 NC=B 
GSURFACE 61 60 70 71 EL1=1 EL2=4 NO=9 NC=B 
GSURFACE 71 70 85 87 EL1=1 EL2=1 NO=9 NC=B 
*NODES FOR THE WELD ARE GENERATED 


LINE 3 100 NCENTER=200 EL=3 M=1 
CINE 100 204 126 152 178 
LINE 204 211 NCENTER=209 EL=1 M=1 


LINE 211 210 NCENTER=209 EL=1 M=1 
Linz 210°39 208207 206 

LINE 
*THE MESH FOR THE WELD IS GENERATED 


A 
N 
A 
LINE C 100 211 204 
A 
N 
C 


GSURFACE 100 3 39 211 EL1=3 EL2=3 NO=9 NC=B 


*NODES FOR THE HEADER ARE GENERATED 
LINE C 100 39 204 2112 210 

LINE S 50 6 EL=3 M=1 

LINE S 6 61 EL=3 M=1l 

LINE C 50 61 6 

*THE MESH FOR THE HEADER IS GENERATED 
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APPENDIX G. ADINA-IN INPUT FILE FOR THE AXISYMMETRIC 


OOp pS 
NUP fb 


oooooojyoamwoo 


NWOORPRPRPNNNNWWN 


am 
a) 
~ 


2297 30eSoo 
2.93059481 
2.88 
2.84814711 
2.844 
2.969 
2.98666666 
2.96 


GSURFACE 50 100 39 61 EL1=2 EL2=6 NO=9 NC=B 

GSURFACE 6 61 72 75 EL1=3 EL2=4 NO=9 NC=N 

*PORTIONS OF THE MODEL ARE DEFINED FOR PLOTTING 

EZONE TUBEO1 

ieorer TO 3 

EZONE TUBE02 

fe tee LO 10 

EZONE WELD 

Tester cierto 19 

EZONE HEAD 

20: Sree 1 TO 43 

BZONE CLOSE XMIN=0 XMAX=0 YMIN=.63 YMAX=1.0 ZMIN=2.782 ZMAX=3.1 
*THE MASTER DEGREES OF FREEDOM ARE DEFINED 

MASTER REACTIONS=YES IDOF=100111 DT=1 NSTEP=2 ITPS56=1 
PRINTOUT MAX 

ANALYSIS TYPE=STATIC 

ATHE LOAD 1S APPLIED IN 2 STEPS 

TIMEFUNCTION 

0 0 

eo 

« 

TOL T=E 1=15 DN=.00001 

*NODES ONLY USED FOR GENERATION ARE FIXED, THE BOTTOM RIGHT HAND 
ANODE OF THE HEADER IS FIXED 

BOUNDARIES 111111 TYPE=NODES /200 209 75 

*THE NODES ON THE RIGHT HAND SIDE OF THE HEADER ARE CONSTRAINED TO 
*ROLL ONLY IN THE Z DIRECTION 

BOUNDARTES 1101ll TYPE=NODES / SO 101 STEP 1 TO 105 6 
BOUNDARIES 110111 TYPE=NODES / 148 156 163 170 177 185 192 
*xTHE CONTACT SURFACE IS DEFINED 

CGROUP 1 CONTACT2 S=AXISYM 

CONTACTSURFACE 1 

Bae 55)162 169 176.1684 191-198 72 

CONTACTSURFACE 2 

G7 56 71 52 48 45 42 38 34 31 61 

CONTACTPAIR 1 T=2 C=l1 F=0 

*THE INTERNAL PRESSURE LOADING 

LOADS ELEMENT 


i. 700 
oleP 1 TO 
roe =1, 700 
fee 700 
372 700 
za 2 700 
ii 2. 700 
Ziez 700 
20 2 700 


*THE LONGITUDINAL TUBE LOAD 
LOADS CONCENTRATED 


Boo 25.5 
S205 =101.3 
BS ga 25.3 


102 346.3666667 
103 344.55 

W045 342.7333333 
LOS 340.916666/ 


ig. 


6 339.8 
148 327.9625 
56 9316-825 
163, 305.6875 
170 294-55 
177 283.4125 
185 272.275 
92) 26.375 


75 250 
a4 350 
(i 350 
10 350 
a3 350 
16 350 
ro 350 
37 350 
23 350 
26 g50 
ZS 250 
60 350 
a3 350 
36 350 
41 250 
bow 350 
47 eo 
ap 350 
54 350 
70 350 
58 350 
85 350 
2 350 
3 250 
62 350 
63 350 
64 350 
65 350 
66 300 
100 350 
23! 7250 
Lie. eoS0 
itl +250 


DATAB CREATE 

*AINPULT FOR THE FINE MESH, 2-DIMENSIONAL AXISYMMETRIC MODEL 
HEAD ‘HEADER’ 

CONTROL B=4800 SYSTEM=2 DEVICE=2 OPTION=1024 

CONTROL PLOTE=WAIT PLOTS=NO 

ACOORDINATES ARE DEFINED 


COORDINATES / ENTRIES NODE xX x vA 

Ak 0 - 63 a0 

3 0 =o70 2.0 
a7 0 63 2.844 
a2 0 S25 2.844 
60 0 oS 227 o2 
61 0 5 Za7eZ 
70 0 Joo 0 
gal 0 SUS 1.0 
V2 0 s 2 O00L 10 
85 0 2os 0.0 
87 0 sf 0.0 
100 0 -938 3.0 
200 0 soe Zena? 
i572 0 Bigs e's 297 soocac 
256 0 ae lelal 2.88 
248 0 ~792 2.844 
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Zod 0 O23 22969 

126 0 .938 2.98666666 
178 0 dae 2.96 

250 0 ol 2.844 

246 0 “i fai a 2-844 

50 0 1.0 oa 

6 0 10 PAP Bo 7 

75 0 1.0 150 

a3 0 2930 2.993395335 
139 0 3956 2298 

165 0 —Jo6 2. 96666666 
249 0 -8025 2-844 

247 0 s7e040 2-844 

245 0 . 7605 2.844 


*AMATERIAL PROPERTIES ARE DEFINED 
MATERIAL 1 ELASTIC E=29600000 NU=.3 
*ATHE ELEMENT IS DEFINED 

EGROUP 1 PLANE SU=AXISYMM 

\ToHe, MESH FOR THE TUBE IS GENERATED 
GSURFACE 3 1 37 39 EL1=2 EL2=6 NO=9 
GSURFACE 39 37 60 61 EL1=2 EL2=4 NO=9 
GSURFACE 61 60 70 71 EL1=2 EL2=8 NO=9 
GSURFACE 71 70 85 87 EL1=2 EL2=2 NO=9 
*NODES FOR THE WELD ARE GENERATED 
LINE A 3 100 NCENTER=200 EL=6 M=1 
ieee N 100 178 113 126 139 152 1695 
LINE A 178 256 NCENTER=251 EL=3 Mel 
mene, © 100 256 178 

LINE A 256 250 NCENTER=251 EL=3 M=l 
LINE N 250 39 249 248 247 246 245 
mene € 256 39 250 

aTHe MESH FOR THE WELD [S GENERATED 
GSURFACE 100 3 39 256 EL1=6 EL2=6 NO=9 
*NODES FOR THE HEADER ARE GENERATED 
mone C 100 39 178 256 250 

LINE S 50 6 EL=6 M=1 

fone S 6 61 EL=6 M=l1 


fone, C 50 61 6 

Sire, MESH FOR THE HEADER IS GENERATED 

GSURFACE 50 100 39 61 EL1=4 EL2=12 NO0=9 

GSURFACE 6 61 72 75 EL1=6 EL2=8 NO=9 NC=N - 

*PORTIONS OF THE TUBE ARE DEF INED-POR~PROTEING- 

EZONE TUBEO1] - 

imocer 1 TO 12 

EZONE TUBE02 

iseoler 1 TO 44 

EZONE WELD 

Soeotrre 1 TO 80 

EZONE HEAD 

geo ltere | TO 256 

BZONE CLOSE XMIN=0 XMAX=0 YMIN=.63 YMAX=1.094 ZMIN=2.688 ZMAX=3.1 
BZONE HEADO] XMIN=0 XMAX=0 YMIN=1.094 YMAX=2.0 ZMIN=2.688 ZMAX=3.0 
BZONE HEADO2 XMIN=0 XMAX=0 YMIN=1.094 YMAX=2.0 ZMIN=1.0 ZMAX=2.688 
BZONE HEADO3 XMIN=0 XMAX=0 YMIN=.63 YMAX=1.094 ZMIN=1.0 ZMAX=2.688 
BZONE TUBE XMIN=0 XMAX=0 YMIN=.63 YMAXK=.693 ZMIN=0 ZMAX=1.0 

@THeE, MASTER DEGREES OF FREEDOM ARE DEFINED 

MASTER REACTIONS=YES IDOF=100111 DT=1 NSTEP=2 

PRINTOUT MAX 

ANALYSIS TYPE=STATIC 

fore LOAD IS APPLIED IN 2 STEPS 


fT MED UNCTION 
0 0 
7. gf 


x 

TOL T=E I=15 DN=.00001 

*ANODES USED FOR GENERATION ARE FIXED, THE BOTIOM RIGHT HAND NODE 
MOF THE HEADER IS FIXED 

BOUNDARIES 111111 TYPE=NODES / 200 251 75 


79 


*THE NODES AT THE RIGHT HAND SIDE OF THE HEADER ARE CONSTRAINED TO 
*ROLL ONLY IN THE Z-DIRECTION 

BOUNDARIES 110111 TYPE=NODES/ 50 368 STEP 1 TO 378 6 558 STEP 13 TO 740 
*THE CONTACT SURFACES ARE DEFINED 

CGROUP 1 CONTACT2 S=AXISYM 

CONTACTSURFACE 1 

87 210 205 199 71 191 186 181 175 170 164 159 154 148 143 137 132 127 T2i)iieees 
CONTACTSURFACE 2 

61 570 583 596 609 622 635 648 661 674 687 700 713 726 739 752 72 
CONTACTPAIR 1 T=2 C=l1 F=0 

*ATHE INTERNAL PRESSURE LOAD 

LOADS ELEMENT 

2 = 700 

Same 62150 


*THE LONGITUDINAL TUBE LOAD 
LOADS CONCENTRATED 

85 3 -12.66666666 

217 ees ~50.66666666 

2 kOe a3 =25e 35553953 

2i5a 43 -50 .66666666 


87 2 -12.66666666 
*THE THERMAL LOAD 
LOADS TEMP TR=0 
50 350 

368 349.0916667 
369 348.1833333 
370 347.275 

37 346.3666667 
372 345.4583333 
373 344.55 

374 343 .6416667 
375 $42.7333333 
376 341.825 

LY 340.9166667 
S73 340.0083333 
6 S39 

258 334255125 
S71 32429025 
584 3226399515 
597 S16-e25 

610 Si. 2oece 
623 305.6875 
636 300.11875 
649 294-00 

662 288.98125 
67> 283.4125 
688 277.84375 
701 PAE PARP AT =, 

714 266. 70625 
727 261.1575 
740 255. 56075 


ID 250 
a 350 
olor i £0 
5 350 
218 350 
See ee co 
zee 350 
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100 250 
390 250 
Siar.) TO 
326 SoG 
50 350 
ae 350 
ster 5 TO 
36 350 
43 350 
48 S50 
54 350 
39 S50 
66 350 
oy 356 
78 Sole 
83 350 
90 350 
95 350 
OL 350 
106 S50 
il 1 350 
60 SoC 
120 a50 
25 350 
31 550 
36 350 
142 S20 
147 350 
153 350 
58 hel 0) 
163 250 
169 350 
174 350 
180 350 
a5 350 
9 0 350 
E35 356 
70 350 
204 350 
Z09 550 
214 Sis\(e) 
85 S50 
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APPENDIX H. ADINA-IN INPUT FILE FOR THE 3-DIMENSIONAL 
CONTACT SURFACE 


ATHESE ARE THE CONTACT SURFACES FOR THE 3-DIMENSIONAL MODEL 
CGROUR 1 CONTSEES 
ACONTACTSURFACE 1 IS THE HEADER AROUND THE INNER TUBE 
COMPACT SUREACE 2 

P25 <bZ4 (970 9963S 

yeg (9/09 982 Boat 

981 “982 “924 255 

2937 (994 69 68 

P24" 125. “27 eee e 


r- 

OWON OM PWN PF 
\O 
~ 
© 
\O 
~ 
ht 
\O 
Co 
GJ 
\O 
OO 
N 
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1169 


she pe) 


ACONTACTSURFACE 2 
CONTACTSURFACE 2 


585 
S23 
oS 1 
339 

586 
1354 
358 
136 2 

587 
SO 
ESS 
ii 96 

588 
Ihikegk 
1189 
ie 7 

589 
rs Z 
1190 
JL IES Rs} 

520 
eZ 2 
1218 
1224 

Bi 1 
23 
E219 
222 

BioZ 
250 
1258 
1266 

S93 
Z> L 
Pasay, 
ZG 7 

594 
1286 
1294 
E302 

473 


26 
roo 
256 
1362 

587 
1180 
IL dbeite: 
1196 

588 
We detdl 
Thats, 
ey 

) eye) 
Tlie 
Maes: 
P38 

590 
IAI 
IL lus) 
1224 

Soe 
eS 
PZ Lg 
I PPA 

Bye) 
Z5 0 
2513 
1266 

393 
Zee! 
Zo 
26% 

594 
1286 
1294 
ies Oz 

473 
1287 
a 7S). 
1303 

584 


1354 
Ieyeys' 
eo Z 
2)8) 1! 
1180 
1188 
ILI ke: 
35.2 
Tiel 
a Liahs, 
Lig? 
533 
Pez 
IGS 19 
iLJUShe 
534 
a eZ 
IPA dks: 
1224 
335 
JANES: 
PALS, 
25 
256 
ZO 
ios 
1266 
5)8) ]/ 
Zou 
IL aeye, 
EZ 
538 
1286 
1294 
iS 02 
559 
1287 
IPSN 
S02 
468 
ILE AZ 


iS THE HEADER AROUND THE OUTER TUBE 


le 2S 
igsieul 
ILS SS, 
250 
1354 
1US oye: 
So 
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E287 
aS ps 
1303 
984 
ig PAA 
£330 
SUS) eye! 


L322 
Too 
36 

Byay 
IEPs) 
bees! 
i359 


T2350 
£338 

a2 
Las 
Yoon 
1339 

5350 


xXCONTACTSURFACE 3 


eo 


WON AW PWN 


NTACTSURPAGE 3 
195 237 238 196 
Za 90S) 9238 
DON 297 72946000 
297 339 340 298 
SIS EN Src Bred 3440 
Breil (ys| 53S) Sis 
196 238 239 197 


255. 9197. oo 
21 2986299 7 


298 
340 
332 


340 341 299 
382 383 341 
675 /O F385 
E97 239240098 
bSo 9s 206 
JL) ZIDeS 00, 9s 
299 341 342 300 


341 383 384 342 
333, /Om7 lL 354 
198 240 241 199 
240 93 94 241 
93 300 301 94 
300 342 343 301 
342 384 385 343 
384 717 2 ast SD 
199 241 242 200 
241 94 95 242 
94 301 302 95 
301 343 344 302 
343 385 386 344 
S85 725735. See 
200 242 243 201 
242 95 96 243 
Son OU Zeo0s 726 
302 344 345 303 
386 387 345 
73 74 387 
243 244 202 
96 97 244 
9653030304 897 
303 345 346 304 
345 387 388 346 
38/7 /4, 75 362 
202 244 245 203 


IS THE INNER TUBE 
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44 244 97 98 245 
45 97 304 305 98 
46 304 346 347 305 
47 346 388 389 347 
48 388 75 76 389 
49 203 245 246 204 
90 245 98 99 246 
eee oo 305 306 99 
92 305 347 348 306 
53 347 389 390 348 
eeeso9 76 77 390 
595 204 246 235 193 
femet6 99 8 235 
eeees? 306 295 8 
per s06 348 337 295 
Boe ot¢8 390 379 337 
mon 390 77 6 379 
meee oS 235 236 194 
Oe 2535 8 89 236 
Beno 295 296 89 
eee2I5 337 338 296 
eemes/ 379 380 338 
Se, 379 6 67 380 
eeel94 236 237 195 
eam2356 89 90 237 
memos 296 297 90 
pe 296 338 339 297 
fpames38 380 381 339 
f7e3930 6/7 68 381 
MeOnlTACTSURFACE 4 IS THE OUTER TUBE 
CONTACTSURFACE 4 
ea 699 700 658 
Be 592 553 7/00 
Boz (99 760 553 
P59 801 802 760 
801 843 844 802 
843 530 531 844 
soo 700 701 659 
BO S53 554 701 
Bes 760 761 554 
feo S802 803 7/61 
802 844 845 803 
12 844 531 532 845 
ioeea9 701 702 660 
mee 701 554 555 702 
meso 4+ 761 7/62 555 
16 761 803 804 762 
17 803 845 846 804 
Beamer. 532 533 846 
meecoO 702 703 661 
Fume OZ 555 556 /03 
Poa 762 763 356 


MmrWOONNAUMN PW Nr 
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APPENDIX I. ADINA-IN INPUT FILE FOR THE 3-DIMENSIONAL 
PROBLEM 


DATAB CREATE 

HEAD ’ 3D-HEADER’ 

CONTROL B=4800 SYSTEM=2 DEVICE=2 OPTION=1024 
CONTROL PLOTE=WAIT PLOTS=NO 

*THE COORDINATES FOR THE INNER TUBE ARE DEFINED 
SYSTEM 1 TYPE=CYLINDRICAL X=9.4 Y=1 221.125 


COORDINATES /ENTRIES NODE R THETA AL 
it os 0 0 
Z Ave) 0 0 
3 pos 0 156 
4 ATs 0 156 
> 63 0 218 
6 Wes: 0 -218 
u 63 0 220 
8 a Tpe 0 2.0 
2 Os 0 4.0 
10 aS 0 4.0 
eal -79001 0O 220 

ioe CYLINDRICAL 1161 

LINE CYLINDRICAL 2261 

iene, CYLINDRICAL 3 3 6 l 

LINE CYLINDRICAL 4461 

five CYLINDRICAL 5 5 6 1 

CINE CYLINDRICAL 6 6 6 l 

ioeNe CYLINDRICAL 7 7 6 1 

fine CYLINDRICAL & 8 6 1 

iene CYLINDRICAL 9 9 6 1 

LINE CYLINDRICAL 10 10 6 1 

LINE CYLINDRICAL 11 11 61 


Boe MASTER DEGREES OF FREEDOM ARE DEFINED 
MASTER REACTIONS=YES IDOF=000111 

PRINTOUT V=MIN I0=0 

ANALYSIS TYPE=STATIC 

VIEW 1 X=1 

VIEW 2 

VIEW 3 
VIEW 4 Y=-1 
VIEW 5 
VIEW 6 Z=-1 

xTHE MATERIAL PROPERTIES ARE DEFINED 
MATERIAL 1 ELASTIC E=29600000 NU=.3 
pide ELEMENT TS DEFINED 

EGROUP 1 SOLID 

BIneostABLE 11234567 8 21 
EDATA/ENTRIES EL TABLE PRINT 


ik i YES 
Sloe alee 
288 1 Sa) 


etme MESH FOR THE 
GVOLUME 10 10 9 7 7 EL1=6 EL2=1 EL3=4 NO=20 
GVOLUME 8 8 7 7 5 EL1=6 EL2=1 EL3=4 NO=20 
SVOLUME 6 6 5 5 3 EL1=6 EL2=1 EL3=1 NO=20 
GVOLUME 4 4 3 3 1 EL1=6 EL2=1 EL3=1 NO=20 
ee COORDINATES FOR THE CUTER TUBE ARE DEFINED 


SYSTEM 2 TYPE=CYLINDRICAL X=9.4 Y=3.0 Z=1.125 


MW U1 oO 
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COORDINATES /ENTRIES NODE R THETA XL 


463 63 0 ~.4564 

464 73 0 =.49 

465 -o3 0 =093 27705 
466 She 0 =, Lileze, 3 
479 -63 180 =. 1036 

490 ais. 180 —207 

467 S65 0 226 

468 o25 y) 2 

469 -63 0 1.64381482 
470 275 0 LPeGlvoa9259 
471 -63 0 4.0 

472 27S 0 4.0 

473 -7 500170 1.61759041 
S83 279001. 180 1. 945157256 
474 .63 30 -.43276688 
475 -63 60 =. 2650 

476 63 90 =a20 

477 -63 Z0 = sLolS 

478 Slee: 150 rh) 2 olen 
480 03 210 Si PAS PAE Nee 
481 -63 240 a Lo Le 

482 oS 270 = PA 

483 -63 300 =n eis) 

484 263 330 -.43276688 
485 SWE) 30 =2.4601865 205 
486 Ages 60 5 Biel 

487 ies 30 pea 

488 Aho. 0) Say =. 

489 eto L350 -.09813465 
491 t= Z0 =, 09813465 
492 Aa Ge) 240 ree = 

493 Se 5) 270 ea 

494 Path 300 =2365 

495 sto 530 -.46186555 
584 -7o001 30 1.63954765 
565 <7 500) 560 1.69353594 
586 ~7a001 90 1.78148148 
S87 -/ S001 ee uco 1.86342702 
588 « 73001)" 4250 1. 92340532 
520 -/50017 220 12692341531 
S91 ~-75001 240 1.86342702 
S92 -72002 270 1. 78148148 
593 «79001 ~ 300 1.69953594 
594 ~75001 330 1.63954765 
496 263 30 ~. 08036908 
497 {65 60 -.04510387 
498 -63 90 -0030693 
499 -63 120 -05124247 
500 Ba so 150 .08650768 
501 -63 180 .09941564 
502 =63 210 - 08650768 
503 ~63 240 .05124247 
504 63 270 0030694 
505 203 300 -.04510387 
506 203 soc =.08036908 
507 fio 30 = 203626202 
508 57S 60 -.05427972 
S09 eh 90 ~0030673 
SLO ap) 120 -06041831 
SLi ar 150 -10240071 
=p hyd BS 180 . L776 732 
S13 Ra ZO . 10240071 
514 eye) 240 .06041831 
S15 Pa. 270 - 0030693 
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516 AY go 
a27 x75 
940 mek, 
541 63 
542 63 
543 -63 
544 63 
545 ~63 
546 63 
547 Bos 
548 -63 
549 63 
550 63 
551 Tis) 
Sy ZS 
523 75 
554 75 
555 f= 
556 ~75 
S57 5 
558 Fae 
S52 ao 
560 75 
562 oo 


LINE NODES 463 463 474 475 

476 STEP 1 TO 484 
LINE NODES 464 464 485 486 

457 Sieewel FO -495 
LINE NODES 465 465 496 497 

AI 8 sree 1) TO 506 
LINE NODES 466 466 507 508 

S09 Sierra TO S17, 
LINE CYLINDRICAL 467 467 6 1 
LINE CYLINDRICAL 468 468 6 1 
LINE NODES 469 469 540 541 

242 sere TO 550 
LINE NODES 470 470 551 552 

S53 StEPw 10-561 
LINE CYLINDRICAL 471 471 6 1 
LINE CYLINDRICAL 472 472 6 1 
LINE NODES 473 473 584 585 

BBbGe ster el. 0 S94 
kTHE MESH FOR THE OUTER TUBE I5 
GVOLUME 472 472 471 471 470 470 
GVOLUME 470 470 469 469 468 468 
GVOLUME 468 468 467 467 466 466 
GVOLUME 466 466 465 465 464 464 
*THE COORDINATES FOR THE HEADER 


SYSTEM N=3 TYPE=CARTESIAN X=9.4 

COORDINATES / ENTRIES NODE X 

xINNER TUBE HEADER NODES 
325 0.0 
926 15 
327 a fae 
928 Z.0 
S2z9 0.0 
930 eye 
S31 5 dl 
932 22.0 


300 -.05427972 
330 = OSG 20 2ie2 
30 166 229867 
60 1. 71264813 
90 1. 78148148 
eZ0 1.85031484 
150 1.9007043 

180 1.91914817 
210 1.9007043 

240 1.85031484 
270 1. 78148148 
300 tn7 264615 
330 1266225867 
eo 1. 63954951 
60 1.69953704 
90 1.78148148 
120 1. 86342593 
150 1.92341346 
180 1.94537037 
Zo 1.92341346 
240 1-86342593 
270 1. 78148148 
300 1.69953704 
830 Leo 3954951 

GENERATED 


469 469 EL1=6 EL2=1 EL3=4 NO=20 
467 467 EL1=6 EL2=1 EL3=4 NO=20 
465 465 EL1=6 EL2=1 EL3=1 NO=20 
463 463 EL1=6 EL2=1 EL3=1 NO=20 
AROUND THE TUBES ARE DEFINED 


Y=0 Z2=0 
Y VA 
O70 0.0 

6 0.0 0.0 

8 0.0 0.0 
0.0 0.0 
O-0 IE LePas, 

6 0.0 ZS 

8 0.0 ele S 
0.0 Paleo 
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933 C.0 0.0 LB Tiga: 
934 56 0.0 2.25 
J35 22.8 0.0 2-20 
936 Z.0 0.0 Tens, 
937 0.0 2.0 0.0 
338 5156 2.0 0.0 
939 ~218 2.0 Ole 
940 2.0 2-9 0.0 
941 0.0 2.0 Hos gape 
942 «Met 2.0 leLZo 
943 2216 2-0 elZS 
944 2.0 2.0 PeL25 
945 0.0 2.0 Piso hea) 
946 256 Zoo Zea 
947 eels 2.0 Zao 
948 Ze 0 2.0 2.209 
*NODES OF THE CENTERS OF THE TUBES FOR MESH GENERATION 
949 0-0 1.0 Leo 
950 ~ 156 i. 0 Dips Pa 
<pepe y2.0 1.0 ree 
Joe Za0 eo toalz> 
953 0.0 3.0 Pel2o 
954 6 3.0 eZ 
355 3 Pale hls, Ie Bl FAs 
356 zZ.0 3.0 elo 
AOQUTER TUBE HEADER NODES 
1230 1.41 4.7 2.25 
L231) 2.206 4.2764467 2.25 
1232 -.171 4.13822335) 2-23 
12335) -. 56 4.0 2-22 
1234 1.41 4.7 Pelz 
T2398 sale 4.2764467 1.125 
1236 -=.171 4.138223358)172> 
12379) -50 4.0 iS 
1238" 1347 4.7 0.0 
12330 4.2764467 0.0 
1240 =.171 4.1382233500—0 
1241 =-.56 4.0 0.0 
*AGENERATION OF THE MESH FOR THE INNER TUBE HEADER 


LINE NODES 123 125 124 
LINE ARC 68 70 NCEN=951 EL=1 M=l NCO=ALL 

LINE ARC 46 48 NCEN=950 EL=l1 M=l NCO=ALL 

LINE ARC 24 26 NCEN=949 EL=1 M=l NCO=ALL 

GVOLUME 948 936 125 123 947 935 70 68 ELl=1 EL2=1 EL3=4 NC=N 
GVOLUME 947 935 70 68 946 934 48 46 ELl=1 EL2=1 EL3=1 NC=A 
GVOLUME 946 934 48 46 945 933 26 24 ELl=1 EL2=1 EL3=1 NC=A 

x 

LINE NODES 125 127 126 

LINE ARC 70 72 NCEN=951 EL=1 M=1 NCO=ALL 

LINE ARC 48 50 NCEN=950 EL=1 M=1 NCO=ALL 

LINE ARC 26 28 NCEN=949 EL=l M=1 NCO=ALL 

GVOLUME 936 932 127 125 935 931 72 70 ELl=l EL2=1 EL3=4 NC=N 
GVOLUME 935 931 72 70 934 93% 50 48 EL1=l EL2=1 EL3=1 NC=A 
GVOLUME 934 930 50 48 933 929 28 26 ELl=l EL2=1 EL3=1 NC=A 

* 

LINE NODES 127 129 128 

LINE ARC 72 74 NCEN=951 EL=1 M=1 NCO=ALL 

LINE ARC 50 52 NCEN=950 EL=l1 M=1 NCO=ALL 

LINE ARC 28 30 NCEN=949 EL=1 M=1 NCO=ALL 

GVOLUME 932 928 129 127 931 927 74 72 ELl=l EL2=1 EL3=4 NC=N 


90 


GVOLUME 931 927 74 72 930 926 52 50 EL1=1 EL2=1 EL3=1 NC=A 
GVOLUME 930 926 52 50 929 925 30 28 EL1=1 EL2=1 EL3=1 NC=A 
re 


LINE NODES 129 131 130 

LINE ARC 74 76 NCEN=951 EL=1l 
LINE ARC 52 54 NCEN=950 EL=1 
LINE ARC 30 
GVOLUME 928 
GVOLUME 927 
GVOLUME 926 
* 

LINE NODES 131 11 132 

LINE ARC 76 6 NCEN=951 EL=1 M=1 NCO=ALL 
LINE ARC 54 4 NCEN=950 EL=1 M=1 NCO=ALL 
LINE ARC 32 2 NCEN=949 EL=1 M=1 NCO=ALL 
GVOLUME 940 
GVOLUME 939 


* 
bane NODES Tr 123 122 

LINE ARC 6 68 NCEN=951 EL=1 M=1 NCO=ALL 
LINE ARC 4 46 NCEN=950 EL=1 M=1 NCO=ALL 
LINE ARC 2 24 NCEN=949 EL=1 M=1 NCO=ALL 


32 NCEN=949 EL=1 M=1 NCO=ALL 
940 131 129 927 939 76 74 EL1=1 EL2=1 EL3=4 NC=N 
939 76 74 926 938 54 52 EL1=1 EL2=1 EL3=1 NC=A 
938 54 52 925 937 32 30 EL1=1 EL2=1 EL3=1 NC=A 


944 11 131 939 943 6 76 ELl=1 EL2=1 EL3=4 NC=N 
943 6 76 938 942 4 54 EL1l=1 EL2=1 EL3=1 NC=A 
GVOLUME 938 942 4 54 937 941 2 32 EL1=1 EL2=1 EL3=1 NC=A 


GVOLUME 944 948 123 11 943 947 68 6 ELl=1 EL2=1 EL3=4 NC=N 


GVOLUME 943 947 68 
GVOLUME 942 946 46 
*GENERATION OF THE 
LINE NODES 587 589 
LINE NODES 532 534 
LINE NODES 510 512 511 

LINE NODES 488 490 489 

GVOLUME 948 944 589 587 
ine NODES 3532 510° 891 

GVOLUME 947 943 534 532 
LINE NODES 510 488 909 

GVOLUME 946 942 512 510 
* 


4 941 945 24 
MESH OR Sine 
588 
323 


947 943 534 532 


946 942 512 510 


945 941 490 488 
mee NODES 589 5919590 

LINE NODES 534 536 535 

LINE NODES 512 514 513 

LINE NODES 490 492 491 

GVOLUME 944 940 591 589 
LINE NODES 536 514 893 

GVOLUME 943 939 536 534 
LINE NODES 514 492 911 

GVOLUME 942 938 514 512 
* 

Give NODES 591 593 592 

LINE NODES 536 538 537 

LINE NODES 514 516 515 

LINE NODES 492 494 493 


$43 939 Ss6 53¢ 


942 938 514 512 


941 937 492 490 


GVOLUME 940 1238 593 591 939 12399538 536 ELIel ELZ=1 EL324 NC=N 


DINE NODES 538 516 894 


EYSLUME 939 1239 538 536° 938912240 516 514 EGl=! ELZ=1 EL3=1 NC=A 


LINE NODES 516 494 912 


6 942 946 46 4 EL1=1 EL2=1 EL3=1 NC=A 
2 EL1=1 EL2=1 EL3=1 NC=A 
OUTER TUBE HEADER 


EL2=1 EL3=4 
EL2=1 EL3=1 
EL2=1 EL3=1 


EL2=1 EL3=4 
EL2=1 EL3=1 
EL2Z=1 EL3=1 


NC=N 
NC=A 


NC=A 


NC=N 
NC=A 


NC=A 


GVOLUME 938 1240 516 514 937 1241 494 492 ELl=1 EL2=1 EL3=1 NC=A 
* 


LINE NODES 593 473 594 
LINE NODES-538 468 539 
LINE NODES 516 466 517 


oh 


. GVOLUME 1231 947 532 


LINE NODES 494 464 495 
GVOLUME 1238 1234 473 593 
GVOLUME 1239 1235 468 538 
GVOLUME 1240 1236 466 516 
x 


LINE NODES 473 585 584 
LINE NODES 468 530 529 
LINE NODES 466 508 507 
LINE NODES 464 486 485 
GVOLUME 1234 1230 585 473 
LINE NODES 530 508 890 
GVOLUME 1235 1231 530 468 
LINE NODES 508 486 908 
GVOLUME 1236 1232 508 466 
x 

LINE NODES 585 95687 7oe6 
LINE NODES 5303522) 531 
LINE NODES 508 s10—502 
LINE NODES 486 488 487 
GVOLUME 1230 948 587 


GVOLUME 1232 946 510 
xCOORDINATES FOR THE 
SYSTEM N=0 

COORDINATES/ ENTRIES NODE 
1366 
1367 
1368 
1369 
roa 
Li3s7s 
1372 
Ls 
1374 
1375 
1376 
WS Ae 
1376 
Nes ys 
1380 
Ser 
sez 
1383 
1384 
i325 
1386 
rSss7 
1388 
oe? 
1390 
jk BE 
1332 
1332 
1394 
1395 
L3o6 
1337 
Tai38 
1339 
1400 


iZs9 
1240 
1241 


L225 
1236 
i227 


468 
466 
464 


iZ35 
£23601 
WHE TE 


508 


486 


585 22510947 5320530 


ys 
GasZ 
O32 
4.67 
4.67 
6.79 
6.79 
4.85 
4.85 
6207 
6.57 
4.6 
4.6 
Sys) Fi2) 
28 2 
4.057 
4.057 
4.61 
4.61 
BIS 
3233 
3.03 


a2 


530 


528 
516 
494 


468 
466 


464 


N 
ul 


NONONONOWN 
e 7 e e ® 6 e 6 
N N) N NJ N N N N 
ul ul uw WwW ul in uw uw 


N 
mn 


e 6 ® ® 
tN N N N N 
163) uw uw ul uw 


OSNONONONONONONONONONONONONO 
nO N 
uw uw 


EL1=1 
EL1=1 
EL1=1 


EL2=1 
EL2=1 
EL2Z=1 


EL2Z=1 
EL2=1 
EL2=1 


EL3=4 
EL3=1 
EL3=1 


EL3=4 
EL3=1 
EL3=1 


NC=N 
NC=A 
NC=A 


EL1l=1 EL2=1 EL3=4 NC=N 
530 1232 946 510 508 EL1=:1 EL2Z=1 EL3=1 NC=A 
508 1233 945 488 486 EL1=1 EL2=1 EL3=1 NC=A 
REST OF THE HEADER ARE DEFINED 


*THE MESH FOR THE 
1220 
1231 
Zoe 
1366 
1370 
1374 
1376 
1382 
1386 
390 


GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 
GVOLUME 


zoe 
i222 
1240 
1367 
oa) 
re > 
1379 
1383 
1387 
1391 


i 
REST 
1232 
VZ232 
Zoo 
1370 
1374 
i378 
1382 
1386 
739.0 
1394 


401 


EZs2 
1240 
1241 
let 
T3755 
179 
1363 
ie yely/ 
1392 
1395 


SB. 06Se4-945 2.25 
OF THE HEADER IS GENERATED 


1366 
398 
1400 
369 
i373 
L377 
1381 
[38s 
i289 
i393 


1367 
Dich bey 
1401 
1368 
i372 
376 
Seo 
1384 
1388 
uSo2Z 


P39 
1401 
1369 
bS7 2 
1376 
1380 
1384 
1388 
T3972 
1296 


Tis96 
1400 
ra6e 
373 
ros 
JESy SUE 
13685 
1369 
1393 
eS 7, 


Ei =2 
EL1=2 
EDl=2 
FL1=2 
EL1=2 
EL1=2 
EL1=2 
EL1=2 
EL1=2 
ELL=2 


EL2=4 
EL2=1 
EL2=1 
EL2=1 
EL2=1 
EL2=1 
EL2=1 
EL2=1 
EL2=1 
EL2=1 


EL3=1 
EL3=1 
EL3=1 
EL3=6 
EL3=6 
EL3=6 
EL3=6 
EL3=6 
EL3=6 
EL3=6 


ARE FIXED 
Jake) Fey: 


CENTER OF THE TUBES 
BOUNDARIES 111111 TYPE=NODES/949 STEP 
care NODE aa IME ORIGIN IS FIXED 
BOUNDARIES 111111 TYPE=NODES/1892 
*THE NODES ABOVE THE ORIGIN ARE ONLY ALLOWED TO ROLL IN THE X-DIRECTION 
BOUNDARIES 011111 TYPE=NODES 
Se2e Lolo 20168 1023 1026 1031 1034 1039 931 1043 930 1046 929 
moses Ogease6 19173 1917 1924 1928 1935 1939 1946 1950 1956 
*THE NODES IN THE Y=0 PLANE ARE CONSTRAINED TO ROLL ONLY IN THAT PLANE 
BOUNDARIES 010111 TYPE=NODES 
eos 04s 926 961 1051 966 1017 1053 1054 973 1059 976 1025 1061 1062 
985 1067 990 1033 1069 1070 997 1075 935 1041 1077 927 1004 1079 934 1044 
POgUmoz6 82009 1082 933 1047 1083 925 1394 1395 1891 1893 16894 1896 1900 
BO tet oOcemendant 905 1907 1911 1912 1924 1915 1916 1908 1922 1923 1925 1926 
uez7 £929 1933 1934 1936 1937 1938 1940 1944 1945 1947 1948 1949 1951 1396 
97 1955 1957 
*THE NODES IN THE Z2=0 PLANE ARE CONSTRAINED TO ROLL ONLY IN THAT PLANE 
BOUNDARIES 001111 TYPE=NODES 
ue Zz 3 + 5 og 
61 83 94 105 Tt 
E63 Go leo el 217 
226 Zee oS 2a i 277 
Zes BeIwe cog) Zoos 334 
343 Bez 5355 sol 397 
409 415 418 421 448 
454 460 463 464 473 
490 Sslz 523 534 604 
610 elo- 6254628 667 
673 688 730 
736 745 787 
yas 805 847 
856 ia 7ol 901 
907 916 
1014 1OZZ 
1045 1063 
1it3e 1142 
1167 P79 
LZ 0 1234 
Zon 1296 
1320 1403 
1445 1454 
1493 1504 
1545 LSs6 
1600 1609 
1652 1662 1666 
17.03 1714 
L733 1766 
1810 Jie ee, 


ATHE NODES AT THE 


50 
iS7 
223 
280 
337 
403 
451 
479 
607 
670 
732 
790 
853 
904 


28 
148 


17 
145 
208 
Za: 
o26 
ood 
442 
471 
598 
655 
Zire ih 
as 
838 
B95 


alt 
142 
Z05 
268 
SAS) 
365 
439 
470 
595 
652 
75 
O72 
825 
B92 


9 10 
139 
199 


265 


7 
ibys 
le? 
250 
307 
370 
427 
466 
556 
637, 
697 
754 
edly! 
880 
941 

1030 

IED PA 

1147 

1tea7 

1Z37 
ie. 

1412 

1463 

rSi> 

S67, 

Tele 

1673 

Zo 

V7 

1828 


8 
ios 
I30 
253 
310 
373 
430 
467 
567 
643 
703 
fies 
ezZ5 
883 
942 

OZ 
1124 
1148 
1190 
PAV =, 
£303 
1415 
1466 
MU sphe, 
P5706 
1622 
1676 
i729 
1780 
Te3z 


6 
iG 
184 
247 
SOI 
367 
424 
465 
545 
631 


1040 
1130 
PES3 
lakers! 
1283 
1307 
1425 
1475 
IS Es, 
Sea 
i633 
1685 
1740 
72 
1843 


1042 
Dise 
1154 
1200 
Zao 
1308 
1429 
1478 
f5o¢ 
1565 
1637 
1688 
1743 
1725 
1847 


2136 
163 
1206 
JEP Ale Be, 
113 
1441 
1489 
1542 
i396 
1648 
1699 
io 2 
1806 
1858 


1134 
LIS9 
1204 
1288 
Je akg 
1437 
1486 
1333 
L592 
1644 
1696 
1749 
1802 
1854 
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1861 1865 1869 1872 1876 1880 1883 1886 1889 1898 1909 1920) 1S93ie942 1222 
*THE INTERNAL PRESSURE LOAD IS APPLIED 
LOADS ELEMENT 


lL =-2e cu 
STEP 1 TO 
120 = 27 00 
55 35700 
SltEP 1 TO 
60 -3 700 
115 2.760 
SteEr fro 
L205 =25 7090 
1264-25700 
SLEEP oO, LO 
LOZ 2235700 
203-2 700 
204 -2 700 
2154-3 7700 
SLE eleeLle 
287 35700 
2166—3 =200 
STEP Seca 
288 =3 700 
x 
CONSTRAINTS 


*THE FOLLOWING CONSTRAINTS ARE NECESSARY BECAUSE COORDINATES WERE 
*DUPLICATED IN THE GENERATION OF THE MESH FOR THE HEADER AT THE 
*BOUNDARY OF THE TUBES 

S66 717 eas 1 

968 2 1019 2 

968 371013943 


190207 loos 1 
102052 20S aiez 
10565) 205 la 
1656" 221031 ¥2 
LOS67sel0S les 
DiZ7 ee Oz. 
LEZ 2 eO9Z2 2 
L127 3 O92 
Lido, 1 22s = 
Lisg, 2 222372 
9675) Tleo ss 

96702 Tis0ez 

96755-11603 

980, 16202721 

S80" 27102 7e2 

SB0.3 2102753 

Le2e- 1 eloes1 
Oza 2 1063.2 
i064 12-1099 1 
1064.2 71099 2 
P64 .5. L099 Ts 
1235.1) 7100 1 
en Gers) 74. lal eyo jer 
pos sc Saha k ee eee 
LiG3ei 1136 1 
i637 2 1Lis6a2 


979 12164 1 
S7 9a 64m 
$793 2164-3 
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292M 11035 71 
992-271035° 2 
99273 1035 3 


bose 11071 1 
HOS6. 2 1071 <2 
10721-1107 1 
LOsZ-2 1107 2 
LO72 3 ligy 3 
pias 1 10s 1. 
2243 2 1108 2 
43 3 1108 3 
i67 1 i144 1 
1167 2 1144 2 
Soi el 1168 1 

So ee, es 2 

Solas elles. 3 

bales: Ea ia es Gy ey aa 
woe 2 Ll Zz 
Bros 3 117773 
we 9 1 1200" 1 
Meg) 2 L200: Z 
Press l 120951 
iZz5.2 1209 2 
i253 1209) 3 
Mette) 1243: 2 
eee. 2a 2 
Pik 3s 1248) 3 
1284 1 1249 1 
1284 2 1249 2 
D264 3 1249 3 
Pe ZOn lt 1285) 1 
me20) 2 1285 2 
Pee Lolszi il 
Meme 2 Loci e 
ieo2 3 132123 
i776 L353 31 
igs 2.1353 2 
Pigea 3s 1353 3 
iazsl 1185 1 
Gz 2 1185 Z 
vez s L185 3 
mie? 1 1216 1 
ie?) 2 L216 2 
mrss) 1 12>: 1 
irs S— 21215) 2 
ilsses L215 3 
Me 1256 *1 
Mz7.2 1256 2 
P27 3s 1256 3 
2925 L257 7) 
ieozec 25792 
2927.5 25473 
ice i293 1 
Haze 2.1293 2 
ise i 13829 1 
weS5o.2) 1329) 2 
is5605 1329) 3 
MiB6e. 1357 1 
tiee 12 1357 2 
MB 3 1357 3 


o5 


MPWNRM WN WYNN BWP Pe 


WNWNWNNWNWNWNWNWNNWNWNUNUNUWNNWNWNWNUNBAWNYPUNPD 


1133 
L193 
Eo 3 
1222 
DAFA 
MFA PAL 
Maral 
J FAGASN 
1264 
1264 
1264 
1265 
1265 
1265 
ey ona 
30. 
1337 
oor 
Naa 7 
361 
a6. 
1361 


OLLOWING CONSTRAINTS DEFINE THE CONTACT SURFACES 


Soe 
Sa A 
S50 
553 
554 
554 


WNRMWNHPNRPUWNRPWNRP WNP Neo PR 


WNWNWNNWNWNWNWNWNNWNWNWNWNWNNWUNWNWNDWY 


96 


ZS 
i250 
1286 
1286 
1287 
Ihe) 9474 
ILE 7574 
Le) Seal 
L331 
i358 
IE ey) 
1188 
1188 
BS 
P89 
HEE Ne, 
2S 
2.6 
IAF al, 
hoo 
1258 
I raeys) 
Laas, 
oo 
1294 
1294 
ie 5 
1330 
i330 
P30 
339 
36Z 
T36Z 
iS 
-1 36 
Lo? 
137 
ros 
1224 
1224 
25 
i220 
1266 
1266 
1267 
i267 
1302 
1302 
i038 
336 
1358 

i23 

23 

124 

124 

2S 

ZS 

iZzG 

iZ6 

rer 


NUYNWNWNWNWNNWNWNWNWNWNNWNWNWNWNWNNWNWNWNWNWNNWNWNWNWNWNNWNYWD 


NUNWNWNWNWNNWNWNWNWNWNNWNWNWNWNWNNWNWNUNWNUNNUWNWNWNWNWNYNWYN WD 
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NWNWNNWNWNWNWNWNNWNWNWNWNWNNWNWNWNWNWNNWNWNWNWNWNNWNWNWNWNWD 
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NWNWNNWNWNNWNWNHNWNWNNWNWNNWNWNNWNWNNWNWNNWNNWNWNWNWNWNNWY WH W 


NWNWNNWNWNNWNYWNNWNWNNWNWNNWNWNNWNWNNWNWNNWNNWNWNWNWNWNNYNWN WY WD 
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720 2a ec 
74625 .1317 "3 
447 e275 2 
Tees 17523 
748 21176 2 
159e 20a 2 
(1496371208 3 
750 2 1245 2 
750 312453 
13879271289 "2 
VORA 375.02 
gcle\ = 5). iUe yc) ee) 
78982 1183-2 
7a9@3 1183 3 
79082 1184-2 
79iece lL Ze OZ 
791° 3 121453 
T3221 253m 
792.0. 25a ee 
S29 3202297 22 
830°2 1333 2 
G30n3 1333435 
G3i 2 Lo ieee 
S3les tigi 
BS252 71132 
83302. 220 e2 
83355 22005 
634-2 126172 
834371261 53 
871.27 1305a2 
G7 202) oad 
872-3 .3413 
G7 352 199 2 
O73 msy 99 53 
874 2 1200 2 
8735 27122602 
8753). 122684 
876 2.126902 
B76 3351269. 3 


APORTIONS OF THE MODEL ARE DEFINED FOR GRAPHICS 
EZONE TUBE1] 
LZ coterie re 60 


EZONE TUBE2 

61 62 STEP 1 TO 120 
EZONE TUBE3 

Zi gtZ2 Stee sl 10 156 
EZONE TUBE4 


LS feiss Stee si tO 1oZ 

ZZONE COMBO] TUBE] TUBE3 

ZZONE COMBO2 TUBE2 TUBE4 

ZZONE COMBO3 TUBE1 TUBE2 TUBE3 TUBE4 

EZONE LOAD 

126 Stee 6.50) 192 203 STEP 12 TO 287. 204 Sier 22510-2240 
2ZONE LOADING TUBE] TUBE2 LOAD 
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